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I. INTRODUCTION 

Since its isolation in 2004, graphene — a single sheet of carbon atoms arranged in a honeycomb lattice — has attracted 
an enormous attention because of its highly unusual electronic and transport properties that are strikingly different 
from those of conventional semiconductor-based two-dimensional electronic systems (for a review see Refs. HH3). It 
has been immediately realized the significance and the potential impact of this new material for electronics. This far, 
it has been demonstrated that the graphene has the highest carrier mobility at room temperature in comparison to 
any known material^. However, graphene is a semimetal with no gap and zero density of states at the Fermi energy. 
This makes it difficult to utilize it in electronic devices such as field effect transistor (FET) requiring a large on/off 
current ratio. The energy gap can be opened in a bilayer graphene by applying a gate voltage between the layers^. 
This gate-induced bandgap was demonstrated by Oostinga et all, and the on/off current ratio of around 100 at room 
temperature for a dual-gate bilayer graphene FET was reported by the IBM£. 

Another way to introduce the gap is to pattern graphene into nanoribbons^^. The conductance of graphene 
nanoribbons (GNRs) with lithographically etched edges indeed revealed the gap in the transport measurement !] 1 1 1 12 . 
This gap has been subsequently understood as the edge-disorder-induced transport gap^— rather than the intrinsic 
energy gap expected in ideal GNRs due to the confinement or electron interactions and edge effects^ During last 
years the great progress has been achieved in fabrication and patterning of the GNRs with ultrasmooth and/or 
atomically controlled edges. This includes e.g. a controlled formation of edges by Joule heating^, unzipping carbon 
nanotubes to form nanoribbonsii, chemical route to produce nanoribbons with ultrasmooth edgesi^ and atomically 
precise bottom-up fabrication of GNRs 19 .. All these advances in nanoribbons fabrication will hopefully enable not 
before long the electronic measurement in near-perfect nanoribbons free from the edge or bulk disorder defects. 

An important insight into electronic properties of graphene and GNRs can be obtained from exact analytical 
approaches. The analytic calculations for the electronic structure of the GNRs have been reported in Refs. [20l — 
Q. The electronic structure of the bilayer graphene was addressed in Refs. I25W29I where the analytical results were 
presented (both exact and perturbative). We are not however aware of analytical treatment of bilayer GNRs (Note 
that a numerical study of the magnetobandstructure of the GNRs was reported in Ref. [13 and [3l|, and the analytical 
and numerical treatment of the edge states in the bi- and N-layer graphene and GNRs was presented in Refs. [29| and 
I32I ). The purpose of the present study is to provide an exact analytical description of the spectrum of 7r electrons in 
the zigzag and armchair bilayer nanoribbons and nanotubes including the dispersion relations, the density of states, 
and the conductance quantization. 

The paper is organized as follows: In order to illustrate our method, in Sec.|TT]we present known analytical results 
for a simpler system, monolayer graphene of the finite size. Subsequently in Sec. IIIII we derive the main analytical 
expressions for the energy spectrum of finite-size structures of bilayer graphene. These expressions are used in 
Sec. HVI to analyze the energy spectrum of various bilayer graphene structures near the Fermi energy. Finally, Sec. tVl 
summarizes our findings. 



II. SINGLE LAYER GRAPHENE 

Analytical expressions for the 7r electron spectrum in GNRs and graphene nanotubes (GNTs), based on tight- 
binding model, were provided in Ref. [23l In this section we will rederive the same expressions in an analytically 
simpler way. Our method more clearly shows the connection between solutions for the infinite sheet of graphene and 
for the finite-size sheet. In addition, simpler method will allow us to derive later on analytical expressions of the 7r 
electron spectrum for more complex systems, bilayer GNRs and GNTs. 



A. Electron spectrum in infinite sheet of graphene 

First we will consider 7r electron spectrum in an infinite sheet of graphene. Hexagonal structure of graphene is 
shown in Fig The structure of the graphene can be viewed as a hexagonal lattice with a basis of two atoms per 
unit cell. The Cartesian components of the lattice vectors ai and a 2 are a(3/2, V3/2) and a(3/2, — \/3/2), respectively. 
Here a « 1.42 A is the carbon-carbon distance^. The three nearest-neighbor vectors are given by Si = 0(1/2,-^/3/2), 
82 = a(l/2, — V3/2), and <5i = a{— 1, 0). The tight-binding Hamiltonian for electrons in graphene has the form 

H Er = -tJ2(4bi + b f j a i ), (1) 
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FIG. 1. (Color online) (a) Honeycomb lattice structure of graphene, made out of two interpenetrating triangular lattices, ai 
and &2 are the lattice unit vectors, and Si, i = 1,2,3 are the nearest-neighbor vectors, (b) Indication of labels of carbon 
atoms in the rectangular unit cell, (c) Brillouin zones for hexagonal unit cell (solid hexagon) and rectangular unit cell (dashed 
rectangle). The Dirac points are indicated by solid circles for the hexagonal unit cell and hollow circles for the rectangular unit 
cell. 



where the operators <ij and bi annihilate an electron on sublattice A at site Rf and on sublattice B at site Rf , 
respectively. The parameter t is the nearest-neighbor hopping energy (t rs 2.8 eV). From now on we will write all 
energies in the units of the hoping integral t, therefore we will set t = 1. Let us label the elementary cells of the 
lattice with two numbers p and q. Then the atoms in the sublattices A and B are positioned at R^ = pa.\ + q&2 and 



Rp q = Si +p&i + qa 2 , respectively. 

The 7r electron wave function satisfies the Schrodinger equation, 

H 1 ^ = . 



(2) 



We search for the eigenvectors of the Hamiltonian (JXJ) in the form of the plane waves (Bloch states) by taking the 
probability amplitudes to find an atom in the sites R^ q and R^ g of the sublattices A and B as 



< 9 = cV k ' R - 



Thus Eq. © yields the eiganvalue equations for the coefficients c A and c B 

- Ec A = c B ^(k) , 
-Ec B = c A 4>(-k) , 



(3) 

(4) 
(5) 



where 

0(k) = e 4k ' 61 + e 4k ' 62 + e 4k ' 63 . (6) 
From Eqs. (|4]) and ([5]) we get the eigenenergies and the corresponding coefficients determining the eigenvectors 

A _ 0(k) 



£(k) = sl l0(k)l, 



£(k)' 



c B = l. 



(7) 



where si = ±1. In the anticipation of the rectangular geometry we introduce dimensionless Cartesian components of 
the wave vector 

K = 3ak x , £ = VSaky (8) 
instead of the wave vector components k x and k y . Then using the coordinates of the vectors 8j we have 



(k) 



and the expression for the eigenenergies becomes' 1 ' 



2e l e cos 



E(k) = siWl + 4 cos 



2,£ 



4 COS [ — | cos 



(I) 



(9) 



(10) 
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For satisfying boundary conditions it is useful to adopt a larger unit cell characterized the same geometry as the 
whole sheet of graphene. Since we are interested in configurations of the graphene with rectangular geometry, we will 
use a rectangular unit cell, as it has been done in Ref. [2J. Such unit cell has four atoms labeled with symbols I , A 
, p , r , as it is shown in Fig. [T}d. The atoms with labels I and p belong to the sublattice A , the atoms with labels 
A and r belong to the sublattice B. The position of the unit cell is indicated with two numbers n and m. The first 
Brillouin zone corresponding to the rectangular unit cell contains the values of the wave vectors k, £ in the intervals 
— 7r < k < 7r , — 7r < £ < 7T. We search for the eigenvectors having the form of plane waves, 

jh = c e ^ m + iKn (ll) 

where a = I, p, A, r. This solution can be obtained from Eq. ([3]) using the equalities 

^ = c B , c p = c A e~ lk Sl , c A = c B e- ,;k - ai , c x = c A e~ l2ak * . (12) 

The Brillouin zones corresponding to hexagonal and rectangular unit cells are shown in Fig. [TJ;. Compared to the 
area of the Brillouin zone of the hexagonal unit cell, the area of the Brillouin zone of the rectangular unit cell is two 
times smaller. Smaller Brillouin zone leads to the appearance of additional dispersion branches. Those dispersion 
branches can be taken into account by using two values of the wave vector k in Eqs. (|10[) and (JTJ), the one with 
— 7r < k < 7r and another obtained replacing k by 2ir + n. Using Eqs. ([7]), (|T2"j) we obtain the coefficients of the 
eigenvectors 

■5=-W-'*#4, « = »«-*»+«>, (14) 

where 

0(K,O = « 3 e- < *+2cos^ (15) 

and S3 = ±1 indicates the dispersion branches that appear due to the smaller Brillouin zone. The equation for the 
energy now becomes 



E(k,0 = 51^/1 + 4cos 2 (0 +s 3 4cos cos (|) . (16) 

This equation has been obtained inS 3 .. Zero energy points of the graphene honeycomb lattice with dispersion relation 
(fTU| are at the points K = (2tt,2tt/3) and K' — (2ir, —2ir/3), where coordinates are given in (k,£) space. K points 
correspond to the corners of the first Brillouin zone. Using the Brillouin zone corresponding to the rectangular unit 
cell, the zero energy points have coordinates (0, and the number of these points is only two, as it is shown in 

Fig. [It. 

Since we will consider finite-size graphene sheets, evanescent solutions become important. Solution exponentially 
decreasing or increasing in the a>direction can be obtained by taking n — i\n\ in Eqs. fj 13[) . (|14|) and (|16[) . whereas 
solution exponentially decreasing or increasing in the y-direction can be obtained by taking £ = The dependency 
of the energy on k when £ = is shown in Fig. [2j We see that the branches with real and imaginary k do not intersect 
at \k\ > 0. 



B. Electron spectrum in various single layer graphene structures 



From the boundary conditions we get restrictions on the possible values of the wave vectors k, £. We will consider 
the structures of graphene that have a set of N rectangular unit cells in the x (armchair) direction and a set of N + 1/2 
rectangular unit cells in the y (zigzag) direction, so that there are AT hexagons along the y axis. Note that rectangular 
unit cell shown in Fig. [TJd extends over the whole hexagon in the y direction, whereas it extends over more that one 
hexagon in x direction. 

Using periodic boundary condition, corresponding to the graphene torus, we get that the possible values of the 
wave vectors k, £ are 



2tt . 
Jf j 
2tt 
N V 



N 
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(17) 
(18) 



5 




FIG. 2. (Color online) Left: dispersion branches of graphene for rectangular unit cell, calculated according to Eq. (|16[) . Right: 
dispersion branches for £ = 0, showing propagating solutions (red solid) and evanescent solutions (green dashed). 



Here [-J denotes the integer part of a number. Thus the spectrum of graphene torus is given by Eq. (|16|) replacing n 
and £ by k v and £y. 

For graphene armchair nanotubes one has the periodic boundary condition in the x direction and the requirement 
"0o,n,r = "00,™,; = ipN+i,n,l — ^7V+i.n,r = for the y direction. Since the energy (|16l) . does not depend on the sign 
of wave vector £, we will search for the eigenvectors of the Hamiltonian ([lj as a superposition of periodic solutions 
Eq. (jnj with £ and 

Wn,a = «c Q (£, K 1/ )e* Tn+i " 1 '" + &c a (-£, ^)e- i « m+iK -™ , (19) 

where k„ is given by Eq. (|18[) and £ needs to be determined. From the boundary conditions we get a system of two 
equations for the coefficients a and b 

ac r ,i(£, k v ) + bCr,i(—£, «„) = , (20) 

oe« Wtv,,^, «*) + &e-^ (AA+1) c r j(-e, *v) = . (21) 

This system of equations has non-zero solutions only when the determinant is zero. From Eqs. (|13j>. (|T4|) it follows 
that the coefficients c r j(£, k) do not depend on the sign of £ and we get the condition sin(£(7V + 1)) = or 

e = ^, (22) 

Additionally, there are two AT- fold degenerate levels corresponding to £ = tt with energies E = ±1. The states of 
those levels have zero wave function amplitudes at the I and r sites. 

For graphene zigzag nanotubes one has the periodic boundary condition in the y direction and the condition 
ipm.o,r = i>m,N+i,l — for the x direction. Similarly as for the armchair nanotubes, the energy (|16j) . does not depend 
on the sign of wave vector k, and we search for the eigenvectors of the Hamiltonian (JT]) as a superposition of periodic 
solutions Eq. (jTTJ) with k and — k, 

*Pm,n, a = ac a (^, K )e^ m+lKn + bc a (^,-K)e^ m - iKn , (23) 

where £j is given by Eq. ()17[) and k needs to be determined. From the boundary conditions we get a system of two 
equations for the coefficients a and b 

ac r (Z j ,K) + bc r fe,-K) = 0, (24) 
ac l {Z jl K)e i *( N+1 ') + ba^-ny-^+V = . (25) 

Using Eqs. (fTB")) . (fT4"f we obtain that non-zero solutions are possible when 

sin(/iiV) / £j 



3in(«(J\r+§)) 



= -s 3 2cos ^ . (26) 
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The possible values of wave vector k should obey this equation. The same condition has been obtained in Ref. [23|. 
Equation (|26[) allows for the imaginary values of wave vector k. The imaginary values appear when £ c < < it and 
S3 = —1, where the critical value £ c = 2 arccos(iV/(2iV + 1)) of the wave vector £ is obtained from Eq. (|26|) setting 
k = 0. In the limit N — > 00 from the condition (l2l)l) with imaginary k and Eq. (fT5| follows that E — 0: edge states 
near zigzag edges in the semi-infinite system have zero energy. 

For N x j\f sheet of graphene open boundary conditions in the y direction are the same as for armchair nanotubes 
and in the x direction are the same as for zigzag nanotubes. Since the resulting conditions for the wave vectors k, £ 
are not coupled, the eigenvector of the Hamiltonian ([T]) is a superposition of four periodic solutions having all possible 
combinations of the signs of k and £ and the possible values of the wave vectors are given by Eqs. (f2"2"j) and (|26p . In 
addition there are two TV-fold degenerate levels corresponding to £ = ir with energies E = ±1. 



III. BILAYER GRAPHENE 

Now we will consider the spectrum of it electrons in bilayer graphene. The tight-binding Hamiltonian for electrons 
in bilayer graphene has the form 

i 

~ l ( a ip b 3>P + & L fli <p) ~ tj - J2( a },l a 3^ + a ],2 a J,l) 1 ( 27 ) 

(i,j),P 3 

where the operators di iP and bi tV annihilate an electron on sublattice A p at site Rj p and on sublattice B p at site 

R,f p , respectively. The index p = 1, 2 numbers the layers in the bilayer system. In the Hamiltonian (|2"T)l we neglected 
the terms corresponding to the hopping between atom B\ and atom B2, with the hopping energy 73, and the terms 
corresponding to the hopping between atom A\ (^2) and and atom B2{B\) with the hopping energy 74. Neglect 
of those hopping terms leads to the minimal model of bilayer graphene^. The parameter t± (t± ~ 0.4 eV) is the 
hopping energy between atom Ai and atom A2 while V is half the shift in the electrochemical potential between the 
two layers. Similarly as for the monolayer graphene, we will express all the energies in the units of t. 



A. Electron spectrum in infinite sheet of bilayer graphene 



We will proceed similarly as in the previous Section and will analyze an infinite system at first. The atoms in 
the sublattices A\ and A 2 are positioned at Rp,q 2 = p&i + qa.2, in the sublattice B\ the atoms are positioned at 
Rpq = Si + pai + q& 2 and in the sublattice B2 the atoms are positioned at = — <5i + pai + q&2- We search for 
the eigenvectors of the Hamiltonian (|27[) in the form of the plane waves. The probability amplitudes to find an atom 
in the sites R P ,g 2 and R P ,g 2 of the sublattices Aj and Bj are 

The coefficients c Ap and c Bp obey the eigenvalue equations 

-Ec Al = Vc Al +c Sl 0(k) + 7 c Aa , (29) 

-Ec Bl = Vc Bl + c Al 4>(-k) , (30) 

-Ec M = -Vc A ' 2 + c B2 4>(-k) + 1C Al , (31) 

-Ec B2 = -Vc B2 + c A2 4>(k) . (32) 

Here energy E, potential V and interaction between layers 7 = t±/t are in the units of the hoping integral t. Using 
the nearest-neighbor hopping energy t w 2.8 eV and the hopping energy between two layers t± sa 0.4 eV one gets 
7 w 0.14. When V = 0, the tt electron spectrum is determined by the equation 



E(k)= Sl [ s '+Jl- + \^ , (33) 
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where st, Sa = ±1. The coefficients of the eigenvector are 



c A 1= _E0± C B 1=1 (34) 

0(-k)' 



A2 g(k) Bq 0(k) 



-k) fl>(-k) 



When V ^ the spectrum is 




r^-|o(kl|- + .s,\/^ + |0(k)P(4^+72) (35) 



and the coefficients of the eigenvector are 



A E(k) + V B , , 

^ = m-v cft = _iw_ /(k) , 

0(-k) v ; ^(-k) 

where the function 

7(J5(k) - V) 

describes the contribution of the second sheet of graphene to the eigenvector. 

Finite-size bilayer graphene sheets can be in AB-q or AB-/3 stacking, as is shown in Fig. [5Ji,b. Similarly as for 
graphene monolayer, we will use rectangular unit cells, one shifted with respect to the other, in each layer of bilayer 
graphene. However, the position of rectangular cells are different for different stacking types. Rectangular unit cells 
have eight atoms with labels It , Ai , pt , r\ and I2 , A 2 , Pi , f% , as is shown in Fig. (3J;,d. For the AB-a stacking the 
atoms with labels It, p\ belong to the sublattice A\ , atoms Ai, T\ to the sublattice Si, atoms I2, pi to the sublattice 
A2 and atoms A2, r2 to the sublattice B2. For the AB-/3 stacking the atoms with labels It, p\ belong to the sublattice 
Bi , atoms Ai, r\ to the sublattice A\, atoms I2, P2 to the sublattice B 2 and atoms A2, r 2 to the sublattice A 2 . 

We search for the solutions of the form 

lpm,n,a p = c a p ^ ' ' (^8) 

where a = I, p, A, r is the label of atoms and p = 1, 2 is the number of the layer. For the AB-a stacking this solution 
can be obtained from Eq. (|2"51) using the equalities 

c ri = c Bl , c Pl = c Al e- lk Sl , c Al = c Sl e- 4k ' ai , c h = c Al e -' l2ak * , (39) 
c T2 = c B2 e~ %ak * , c P2 = c A2 e- %k Sl , c A2 = e B V k ' 52 , c h = c M e lak * (40) 

whereas for the AB-/3 stacking the coefficients are 

c n = (c Al ) m , c pi = (c B ^)*e- ik ^, c Xl = (c^)*e- <k - ai , c h = (c B ^)* e~ i2ak % (41) 
c r2 = (c A =)*e- ik ' a = , c P2 = (c S2 )*e-^, c A2 = (c^)* , q a = (c s =)V k ^ . (42) 

Similarly as for monolayer graphene, to take into account the smaller Brillouin zone we need two dispersion branches: 
one with k and one with 2n + k . Using Eq. (jM)) or Eq. (|3"r?)) we obtain the coefficients of the eigenvectors. The 
expressions for the coefficients are presented in Appendix [A] The expression for the energy becomes 



E(k, t)= s Jf- + V 2 + \cj>( K , Ol 2 + s 2 y ^ + |0(«, 01 W 2 + 7 2 ) (43) 



which reduces to 



7 , /7 2 



B(K,e) = Si S2 ^ + ^ + |0(k,OI 2 (44) 
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FIG. 3. (Color online) Upper part: sublatices A% , A% , Bi , B2 on bilayer graphene in AB-a stacking (a) and AB-/3 stacking 
(b). Lower part: indication of labels of carbon atoms used in the description of the ir electron spectrum for the bilayer graphene 
with AB-q stacking (c) and AB-/3 stacking (d). 



for V = 0. Here 

|0( K ,O| 2 = l + 4cos 2 (0 +S3 4cos(0cos(|) (45) 

and S3 = ±1 indicates the dispersion branches that appear due to the smaller Brillouin zone. 

In addition to the propagating waves, for finite-size bilayer graphene sheets evanescent solutions become important. 
Solution exponentially decreasing or increasing in the x-direction can be obtained by taking k = i\n\. Solution 
exponentially decreasing or increasing in the y-direction can be obtained by taking £ = In addition to the purely 
imaginary £ there are solutions, corresponding to S3 = — 1, having complex values of £. The dependency of the energy 
on the wave vector n when the wave vector £ is constant and on the wave vector £ when the wave vector k is constant 
is shown in Fig. 2] We see that now, in contrast to the graphene monolayer, the branches with real and imaginary k 
can have the same energy. 



B. Electron spectrum in various bilayer graphene structures 



We will consider the structures of bilayer graphene that have a set of N rectangular unit cells in the x (armchair) 
direction and a set of Af + 1/2 rectangular unit cells in the y (zigzag) direction, so that there are AT hexagons along 
the y axis. Note that the rectangular unit cells shown in Figs. [3t and EJi extend over the whole hexagon in the y 
direction, whereas they extend over more that one hexagon in x direction. In principle, in the case of bilayer graphene 
nanotubes the numbers N or Af in for the inner and outer cylinders arc different. However, for simplicity we will 
consider them as the same, which is a good approximation for sufficiently large tubes when N — > 00 or Af — > 00. 

Similarly as for graphene monolayer, from the boundary conditions we get restrictions on the possible values of the 
wave vectors k, £. Using periodic boundary condition, corresponding to the bilayer graphene torus, we get that the 
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FIG. 4. (Color online) Dispersion branches of bilayer graphene: dependency of the energy on the wave vector re when the wave 
vector £ is constant (£ = 0) (left) and on the wave vector £ when the wave vector re is constant (re = 1.0) (right). Propagating 
solutions are shown with red solid line, evanescent solutions with green dashed line, and evanescent oscillating solutions with 
complex value of the wave vector £ are shown with blue dotted line. In order to show the structure of the dispersion branches 
more clearly, the value of the parameter 7 is set sufficiently large, 7 = 0.5. 



possible values of the wave vectors k, £ are given by Eqs. (TT7|) . (fl"5j) . 

For bilayer graphene armchair nanotubes one has the periodic boundary condition in the x direction and the 
condition 

V'0,n,r p = ^0,n,i„ = "VV+l.rUp = VW+l,n,r p = (46) 

for the y direction. Here p = 1, 2 is the number of the layer. This condition is the same for both the AB-a and AB-/3 
stackings. For bilayer graphene with AB-a stacking the coefficients c rpj ; p (£, re) do not depend on the sign of £ and we 
get the same conditions (fT5|) , (|2"2"|) for the wave vectors k, £, as for the monolayer graphene armchair tubes. 

For bilayer graphene with AB-/3 stacking the coefficients c rpi ; p (£, k) depend on the sign of £, and condition for the 
possible values of the wave vector £ is much more complicated. There are eight boundary conditions in the y direction. 
In bilayer graphene there are four eigenstates with different wave vectors along y direction, £W, £( 3 ) and 
having the same energy: i?(re,£W) — ^(k,^ 2 ') = _E(k,^ 3 ') = £?(k, £( 4 '), as is evident from Fig. 0] Two or four of 
the wave vectors fW, £ (2) , f( 3 ), £ (4) can be imaginary or complex numbers. Since the energy does not depend on 
the sign of £, we can form a wave function from superposition of eight waves. From the boundary conditions (|46p 
resulting resulting set of linear equations can have nonzero solution only if 8 x 8 determinant is zero. Analytical form 
of this condition in is too large and too complicated to be useful. 

For bilayer graphene zigzag nanotubes one has the periodic boundary condition in the y direction and the condition 

1pm,0,n = lpm,N+l,li = i J m,0,h = lpm,N+l,r 2 = (47) 

for the x direction. Here p — 1, 2 is the number of the layer. This condition is the same for both the AB-a and AB-/3 
stackings. In the bilayer graphene there are two eigenstates with wave vectors along x direction, k*- 1 ' and t6 2 \ having 
different absolute values but corresponding the same energy: E{k^\ £) = E(k^ 2 \^). One or both of the wave vectors 
k*- 1 ' , re( 2 ) can be imaginary. The energy can be equal only if the signs s±, S2 obey the condition 

4 2) 4 2) = - W (48) 

When the bias potential is zero, V = 0, from the equality of the energy we can express re' 2 ); 

V 2) \ m /«« 



When V ^ then 



2 cos 



(I) 



E(kW,0 (49) 



, 3 cos( — ) = s ^cos( — )± <^AE 2 V 2 +1 2 {E 2 -V 2 ) (50) 

' 2 cos ( 2" 



There are four boundary conditions in the x direction. Since the energy does not depend on the sign of k, we can 
form a wave function from superposition of four waves. From the boundary conditions (|47[) resulting set of linear 
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equations can have nonzero solution only if 4 x 4 determinant is zero. The possible values of the wave vector £ is 
given by Eq. (|17[) , and the conditions for the possible values of the wave vector k are given in the Appendix iBl 

For N x AT sheet of bilayer graphene open boundary conditions in the y direction are the same as for armchair 
nanotubes, Eq. (|4l)|) and in the x direction are the same as for zigzag nanotubes, Eq. (PfTl) . For AB-a stacking, the 
conditions for the possible values of the wave vectors n, £ are combination of the conditions for zigzag and armchair 
bilayer graphene tubes. Specifically, when V = 0, the conditions are given by Eqs. (|22|) and (|B1[) or (|B2|) . When V ^ 
then the conditions are given by Eqs. (|22|) and (|B3|) . For AB-/3 stacking it is impossible to separate conditions for the 
wave vector £ from the conditions for the wave vector k. The resulting expressions are very large and complicated. 



C. Summary of the possible values of wave vectors 

For structures of bilayer graphene, the energy spectrum is completely determined by Eq. (|44|) or (J43J) with appro- 
priate expressions for wave vectors n and £. 

Equations presented in Appendix [B] make one quantum number dependent on the other. This dependence appears 
because of zigzag-shaped edges. For structures where zigzag edges do not exist or their effect can be disregarded the 
wave vector n v can be replaced by a continuous variable. 

Thus, the possible values of wave vectors for various structures are as follows: 

• For the armchair bilayer graphene ribbon of infinite length with AB-a stacking, the wave vectors are determined 
by 



< K < 7T , £j 



(51) 



• For the armchair bilayer graphene ribbon of infinite length with AB-/3 staking we have < n < n and the 
equation for the possible values of £ is complicated. 

• For the zigzag bilayer graphene ribbon of infinite length we have < £ < 7r, the conditions for the possible 
values of k, given in Appendix [B] are different for AB-a and AB-/3 stackings. 

• For the zigzag bilayer carbon tube of infinite length with AB-a or AB-/3 stacking, the wave vectors are determined 
by 



< K < 7T, 



t - — ■ 



3 







N 


+ 1,..., 


N-l 


2 _ 




_ 2 _ 


2 



(52) 



• For the armchair bilayer carbon tube of infinite length with AB-a or AB-/3 stacking: 

2n 



N 




N 


+ 1,..., 


N- 1 


~2_ 




~2_ 


2 



< £ < 7T . 



(53) 



Taking into account the ranges of the possible values of the wave vectors, zero energy points for various structures 
with bias potential V = are as follows: 

• For zigzag bilayer carbon tube zero energy points are (0, =^-j, (0, — ^p). 

• For armchair bilayer carbon tube, zero energy point is (0, ^J. 

• The dispersion of armchair bilayer graphene ribbon has only one zero-energy point (0, 4p) 

• For zigzag bilayer graphene ribbon dispersion this point cannot be shown in the real plane. 



IV. BAND STRUCTURE NEAR THE FERMI ENERGY 



In this Section only a part of the spectrum with smallest absolute value of the energy is in focus. This part 
corresponds to S2 = — 1, S3 = — 1. In order to obtain an approximate expression for the energy spectrum near the 
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Fermi energy we expand Eq. (|45l in power series near the zero point n = 0, £ = 2tt/3, yielding 



V3 



2-y/3 



(54) 




Here q = £ — 2n/3 and \n\ <C 1, \q\ <C 1. Substituting Eq. ([5"4"]) into Eq. (gj) or Eq. (|4"5)) one obtains the approximate 
expression for the energy spectrum. Thus, when the bias potential is zero V = 0, the approximate expression for the 
energy is 



5) = «1 ( ~ 



\ 



7^ 6 / 

T + 7 T+l! 




Furthermore, assuming that |^(k,^)| 2 <C 7, the branch of Eq. (|4"4"]) with S2 = — 1, S3 = — 1 takes the form 



When V 7^ and V < 7, then Eq. (j43|) becomes 



IV 

siV~ Sl — Mk,o\ 2 

7 



si- 



| 0(«,Q| 4 
2F7 2 



(55) 



(56) 



(57) 



The bilayer graphene has a gap at |</>(k, £)| 2 = 2V 2 . However, since the parameter 7 is small, 7 <C 1, the approximate 
expressions (j56")) . ([57]) are suitable only for very small values of |«| and \q\. 

The spectrum of various structures of bilayer graphene can be obtained from the approximate expressions for 
\4>(k, £)| 2 near zero points. In contrast to Eqs. (|56| and (|57|) . the energy of monolayer graphene is E(k,£) = 
s i\f\4>{ K i Ol 2 - Thus, the analysis of the square root of |0(k,£)| 2 essentially was done in Ref. Going back to 
the original wave vectors k x and the band structure of bilayer graphene tubes and ribbons when V = similarly 
as in Ref. [23| can be summarized by the equation 



E v {h\ 



si 



2 
"2 



9 



a 2 [(kn 



(58) 



where fey and k± u denote the longitudinal (continuous) and the transverse (quantized) components of the wave vector, 
respectively. Index a specifies the structure. Further in this Section we will consider only the case when V = 0. 



A. Quantum conductance 

Within the framework of the Landauer approach^—, the zero-temperature conductance of a ideal wire is equal to 

G(E) = ^Y,9»TAE) (59) 

u 

where 2e 2 /h is conductance quantum, <?„ is the band degeneracy, and transmission coefficient T v is zero or unity 
depending on whether the u-th band is open or closed for charge carriers with energy E. 

When bias potential is zero, V = 0, the transmission coefficient is T V (E) = Q(E — E°) for conduction bands and 
T„{E) — Q{\E — E°\) for valence bands. Here E° are the subband threshold energies and Q(x) is the Heaviside step 
function. When the approximation Eq. (|58|) is valid, the subband threshold energies are 

K - s 1 U + ^t + ^kf^ . (60) 

The degeneracies are shown in Table HI The values of g v for armchair bilayer carbon tube and zigzag bilayer graphene 
ribbon with v > 1, represented in Table fl] should be doubled, because electron or hole states with ±fcn 7^ are 
degenerate. 
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TABLE I. Degeneracy gZ of the z/-th band energy \EZ{k\\ = 0)| 



(7 


gZ 


armchair bilayer carbon tube 


= 0), 2(1/^0) 


zigzag bilayer carbon tube 


2 


armchair bilayer graphene ribbon 


1 


zigzag bilayer graphene ribbon 


1 



The electron or hole conductance of armchair and zigzag bilayer carbon tubes and their parent graphene ribbons has 
thus the form of a ladder, symmetrically ascending with the increase in energy for electrons, and with the decrease of 
energy for holes. For the charge carrier energy that falls between the n-th and (n + l)-th bands, the wire conductance 
equals 

armchair bilayer ribbon 

zigzag bilayer ribbon 

. (61) 
zigzag bilayer carbon tube 

armchair bilayer carbon tube 

The conductance for bilayer graphene ribbons has been numerically calculated in Ref. H(J The expression (|6lj) for the 
conductance coincides with that of Ref. H3- 



G{E) 




B. Density of states 



The density of states (DOS) of a quantum wire, including a factor 2 for the spin degeneracy, reads 

v v II ' 

The summation includes all transverse modes with energy E v < E. Using Eq. ([55)1 we obtain the DOS of bilayer 
graphene 

p{E) = J_Y gij ffl+rj e(\E\-\K\), (63) 

3 ™ u va^i - mm\ + m + 7) 

The index v is v — 0, ±1, ±2, . . . for bilayer graphene tubes and armchair bilayer graphene ribbons and v = 0, 1, 2, . . . 
for the zigzag bilayer graphene ribbons. The electron density at zero temperature is obtained by the integration of 
the DOS from the charge neutrality point /io = to the Fermi energy, 

E F 

p{E)dE (64) 
Using Eq. (|6U|) we get 

n a (E F ) = ^-Y.9»V(\Ef\ - \E°M\E F \ + \E°\+t)Q(\E f \ - \E V \) (65) 



C. Armchair bilayer carbon tube 

For the armchair bilayer carbon tube we have that k in Eq. ([54]) has discrete values n u — 2tti//N, v = 0, ±1, . 
and q is continuous. Thus the energy spectrum has the form 



E v {k y ) = si 



2 

2 



( ky ky jj 



2 4ir 2 v 2 
' + 9a 2 A^ 2 



(66) 
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FIG. 5. (Color online) Band structure of armchair bilayer carbon tubes (left), DOS (center) and conductance (right). The 
number of rectangular unit cells in the x direction N = 100. Solid red lines are calculated according Eqs. (|44l) . (|45[) ; dashed 
green lines represent approximation (|66p . Bands with S2 = +1 are not shown. DOS is given in units of a -1 and is calculated 
according to Eqs. (|63[) . ()68[) . The conductance is given in units of 2e 2 /h and is calculated using Eq. (|61l) . 



with 



-n-2,,2 

7T ^ 



Conduction (valence) band bottoms (tops) are equal to 



^ = Si ( v -2 + Vt + ^J (68) 

The distance between the minima of the dispersion branches with the indices s 2 = — 1 and s 2 = +1 is 7. The number 
v of subbands with the index s 2 = — 1 and energy smaller than the energies of the subbands with index s 2 = +1 is 
greater than 1 only when the number of rectangular unit cells in the x direction TV is sufficiently large. Using the 
equation -E„, max = 7 and estimating the subband threshold energy (|68|) as E v « ir 2 v 2 /(jN 2 ) one obtains that the 
requirement !/>l leads to N ^> ir/j. In calculations we used N = 100. 

The band structure, calculated with the use of exact Eqs. (|44|) . (|45|) and approximated according to Eq. (|66| . is 
represented in Fig. [5] One sees that Eq. (|66|) provides accurate reproduction of exact results. Also is shown the DOS, 
calculated with Eqs. (fB5J). ([55]), and the conductance G{E) = (2e 2 /h)2(2n + 1). 



D. Zigzag bilayer carbon tube 



Distinct from armchair bilayer carbon tubes, which are always metallic when V — 0, zigzag bilayer carbon tube has 
a gapless spectrum if j* = A/"/3 is an integer, Ej—j* [k = 0) = 0. Otherwise, zigzag bilayer carbon tube spectrum has 
a gap. If A/"/3 is not an integer, band index of the lowest conduction (highest valence) band can be equal either to 
j* = (Af— 1)/3 or to j* = (VV+l)/3. As a result of expansion near zero-energy points in powers of k and 2ir(j — j*)/Af, 
we arrive at 



K 

"3 



(69) 



where 



j v — || <C 1 , semiconducting 
2 ^ d ( 1 + 5^)« 1 » metallic 



with ^ = 0, ±1, .... The wave vector component n is continuous. The energy spectrum has the form 




3L 

3a 2 



(70) 
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FIG. 6. (Color online) Upper part: band structure of metallic zigzag bilayer carbon tubes (left), DOS (center) and conductance 
(right). The number of hexagons in the y direction M = 102. Lower part: band structure of semiconducting zigzag bilayer 
carbon tubes (left), DOS (center) and conductance (right). The number of hexagons in the y direction M = 100. Solid red 
lines are calculated according Eqs. (1 4411 . (|45[) : dashed green lines represent approximation (|70[) . Bands with S2 = +1 are not 
shown. DOS is given in units of a -1 and is calculated according to Eqs. (|63[) . (|71[) . The conductance is given in units of 2e 2 /h 
and is calculated using Eq. (|61|) . 



Conduction (valence) band bottoms (tops) are equal to 




The number v of subbands with the index s 2 = — 1 and energy smaller than the energies of the subbands with index 
Si = +1 is greater than 1 only when the number of hexagons in the y direction Af is sufficiently large. Approximating 
Eq. dTTJ as 37r 2 t/ 2 /(7VV 2 ) we get Af > v^tt/t- In calculations we used Af = 102 for metallic tubes and Af = 100 for 
semiconducting tubes. 

The band structure, calculated with the use of exact Eqs. (|33|). (|4"5|) and approximated according to Eq. (fTOf for 
metallic and for semiconducting tubes is shown in Fig. [6] One sees that Eq. (1701) provides accurate reproduction of 
exact results. Also is shown the DOS, calculated with Eqs. ([55]) . (fTTj) . and the conductance G(E) = (2e 2 /h)2n. 



E. Armchair bilayer graphene ribbon 



For the armchair bilayer graphene ribbon with AB-ev stacking, the condition for the wave-vector component £ has 
a simple expression. When V — and j* = 2(Af + l)/3 is an integer then the armchair bilayer graphene ribbon is 
metallic. Then index v = j — j* =0 corresponds to the zero-energy band. If 2(Af + l)/3 is not an integer, armchair 
bilayer graphene ribbon spectrum has a gap, and the band closest to zero is cither j* = (2Af+ 1)/3 or j* = (2A/"+3)/3 
depending on which of these two numbers is an integer. For k, v /Af <C 1 we get Eq. (|69[) with 

aTTT I v ~~ \ I i semiconducting 

^( 1 + i7#TT))« 1 ' metallic (?2) 

with v = 0, ±1, .... The wave vector component k is continuous. The difference between the boundary conditions for 
armchair bilayer graphene ribbons and zigzag bilayer graphene tubes results in about two-times smaller band spacing 
in the armchair ribbon spectrum than it was found for the zigzag tube spectrum. 
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FIG. 7. (Color online) Upper part: band structure of metallic armchair bilayer graphene ribbon with AB-a stacking (left), 
DOS (center) and conductance (right). The number of hexagons in the y direction M = 101. Lower part: band structure 
of semiconducting armchair bilayer graphene ribbon with AB-a stacking (left), DOS (center) and conductance (right). The 
number of hexagons in the y direction M = 100. Solid red lines are calculated according Eqs. (I44|) . (145[l : dashed green lines 
represent approximation (|73[) . Bands with S2 = +1 are not shown. DOS is given in units of a -1 and is calculated according to 
Eqs. (|63[l . (|74p . The conductance is given in units of 2e 2 /h and is calculated using Eq. (|6ip . 



The energy spectrum has the form 



Conduction (valence) band bottoms (tops) are equal to 




E„ 



(73) 



(74) 



The number v of subbands with the index = — 1 and energy smaller than the energies of the subbands with index 
Si = +1 is greater than 1 only when the number of hexagons in the y direction J\f is sufficiently large. Approximating 
Eq. ((711) as 37T 2 z/ 2 /(47(A/'+ l) 2 ) we get JV > V3tt/(2j). In calculations we used Af = 101 for metallic ribbons and 
N = 100 for semiconducting ribbons. 

The band structure, calculated with the use of exact Eqs. (|44|) . (|45|) and approximated according to Eq. (|73|) for 
metallic and semiconducting ribbons is shown in Fig. [71 One sees that Eq. (|73p provides accurate reproduction of 
exact results. Also is shown the DOS, calculated with Eqs. (|63|) . (|74"J) . and the conductance G(E) — (2e 2 /h)n. 

For the armchair bilayer graphene ribbon with AB-/3 stacking there are no explicit expressions for the possible 
values of q. 



F. Zigzag bilayer graphene ribbon 



In zigzag bilayer graphene ribbons the wave vector component q is continuous while the possible values of n are 
given by the solutions of the equations (|4"5)) and (|B1I) . (|B2[) for AB-a stacking or Eq. (IB5[) for AB-/3 stacking, presented 
in Appendix [Bj Since equations for k depend on the value of the wave vector £, in zigzag bilayer graphene ribbons 
longitudinal and transverse motions are not separable. When the energy of the subband with the index S2 = — 1 is 
smaller than the energies of the subbands with index S2 — +1, only one of wave vectors and k*- 2 ' is real valued. 
Therefore, the energy subbands can be labeled by the value of k' 1 ' = n v only. 
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FIG. 8. (Color online) Band structure of zigzag bilayer graphene ribbon with AB-a stacking (left), DOS (center) and con- 
ductance (right). The number of rectangular unit cells in the x direction N = 60. Solid red lines are calculated according 
Eqs. (144 1) . Q45p with the allowed values of the wave vector k obtained solving Eqs. (|B1|) . (|B2[) and (|49p ; dashed green lines 
represent approximation (|75|l . Bands with S2 = +1 are not shown. DOS is given in units of a~ . The conductance is given in 
units of 2e 2 jh and is calculated using Eq. (|61[) . 



Depending on the value of |£|, wave vectors k v with v — 0, 1 can become imaginary. There are two critical values 
£ C (i) ; £ C (2) (£ C (i) < £ C (2)) of the waye vector ^ obtained by solving Eqs. (03) and (jBT]), fB2]) or fB5]) with K W = 0. 
When £ C W < |£| < £ c ( 2 ) then one solution kq becomes imaginary whereas in the case £ c ( 2 ) < |£| two solutions Ko and 
Ki become imaginary. Both critical values obey the inequality > 2n/3 and tend to the limit 2tt/3 as the 

number N grows. More tight lower bound of critical values is ^ c ( 1 ' 2 ) > 2 arccos[A^/(2iV + 1)]. The imaginary solutions 
k v represent edge states in zigzag bilayer graphene ribbons. 

The energy bands of bilayer graphene are asymmetric near the point q = 0. The subbands (except corresponding 
to edge states) can be approximated by taking k„ irvjN and the minimum of of the subband located at ^ 2 >: 



Eu{k y ) = S X i-^ + d^r 




(ky ky) 2 + 



9a 2 N 2 



with 



c(2) 



V3a 



(75) 



(76) 



The critical value £ c ^ of the wave vector tends to the limit 27r/3 as the number N grows. Conduction (valence) band 
bottoms (tops) are equal to 



E u 



2 



AN 2 



(77) 



The number v of of subbands with the index S2 — — 1 and energy smaller than the energies of the subbands with 
index S2 — +1 is greater than 1 only when the number of rectangular unit cells in the x direction N is sufficiently 
large. Approximating Eq. ([77|) as ir 2 v 2 /(AjN 2 ) we get N ^ 71-/(27). I n calculations we used N = 60. 

The band structure for zigzag bilayer graphene ribbons with AB-a stacking, calculated with the use of exact 
Eqs. (|44l) . (|45|) with the allowed values of the wave vector k obtained solving Eqs. (|B1|) . (|B2I) and (|49l) . as well as 
approximation (jT5"j) are represented in Fig. [5J Also is shown the DOS and the conductance G(E) = (2e 2 //i)(2n + 2). 
The DOS is calculated from exact band structure and also using Eqs. (|63[) . (|74|) . taking the threshold energies for 
v = 0, 1 to be E v= q \ = 0. The band structure of zigzag bilayer graphene ribbons with AB-/3 stacking is very similar 
to the band structure of ribbons with AB-a stacking, only the critical values £, ctyl \ £, c ^ are slightly different. 

Taking_the limit N — > 00 in the zigzag bilayer graphene ribbon with AB-/3 stacking one can obtain the edge states 
of Ref. [23. However, care should be taken not to loose any solutions. For large number N we can write the absolute 
value of the imaginary wave vector as \k\ = + 5, where is the solution of the equation 



= 2cos(£/2) 

and (5 is a small correction. From Eq. (j4"4"]) it follows that such a value of ensures the equality E(iK^°',(,) 
Expanding Eq. (|44|) in powers of S we get the approximate expression for the energy 



E 



s\5 

1 



2cos 2 f^ -- 



(78) 
= 0. 

(79) 
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There are two eigenstates with wave vectors and having different absolute values but corresponding the same 
energy. From Eqs. (|75|) and (|48p it follows that the corrections to the wave vector obey the condition 

5^ = -ad) . (80) 

In Eq. (|B5I) taking into account only the first-order terms with respect to 5 one obtains the value of the correction 

S = ±2e- 2 ^ N (l~e-^'). (81) 

This expression for the correction is the same as for the single sheet of graphene. The correction 6 decreases expo- 
nentially with increasing the number N. 

For large N it is sufficient to form the wave function obeying boundary conditions (|47p as a superposition of two 
exponentially decreasing terms with the wave vectors and i|«;^ 2 ^|, 

=a< 1 >c^i|K< 1 >^^ (82) 

Substituting this expression for the wave function in the boundary conditions, using Eqs. (|A9[) (|A12|) and taking the 
limit N — > oo we obtain two solutions for the coefficients a^: = and = —a^\ The wave function 
corresponding to the solution a' 2 ' = is localized on the first layer, with the nonzero coefficients and c pi . The 

wave function corresponding to the solution = —a^ contains the difference e~' K< — e~' K< Expanding to 
the first order of 8 we get 

e-I^V _ e -l« <2, l» = e ~(^»-8)n _ e -(«C"»+«}„ w 28ne^ { ° )n . (83) 

Taking the limit N — ¥ oo and dropping the coefficients of the wave function that are of the order of 5 we obtain that 
nonzero coefficients are ip mn ,h, ipm,n,pi m the first layer and ipm,n,r 2 i VVn,n,A2 m the second layer. The coefficients 
i > m,n,r 2 i ipm,n,\ 2 are proportional to e~ K< >n while the coefficients ipm,n,h, ipm,n,pi have ne _K< >n behavior, as in Ref . |29|. 



V. CONCLUSIONS 



An exact analytical description of it electron spectrum based on tight-binding model of bilayer graphene has been 
presented. The bilayer graphene structures considered in this article have rectangular geometry and finite size in 
one or both directions with armchair- and zigzag-shaped edges. This includes bilayer graphene nanoribbons and 
nanotubes. The exact solution of the Schrodinger problem, the spectrum and wave functions, has been obtained and 
used to analyze the density of states and the conductance quantization. Our method brings a connection between it 
electron spectrum in infinite and finite-size bilayer graphene. 
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Appendix A: Eigenvectors of bilayer graphene using rectangular unit cells 



The expressions for the coefficients of the eigenvectors are: For AB-a stacking, V = 

Cri =l, c pi =-e~^§^- (Al) 

c h = -„e-*«-^0 c Xl = S3 e-^"+« , (A2) 

c r2 =-s 1 s 2 —, -r, c P2 = Sl s 2 e 2— -r, A3) 

0(-/C,O 

Cfa =sis 2 s 3 e l * v ' s \ , c A2 = -s lS2 s 3 e l ^ K « / s ; . A4 



For AB-a stacking, V ^ 

, S E + V 



c ri — 1 5 Cp 1 — e ^ 



E-\-V _,'I^ K _]_i1 
0(-K,£) 

Cr 2 = "77 tt/^'O) c P2 =sis 2 e 2 -- / Kj £ ) ) 



For AB-/3 stacking, F = 



For AB-/3 stacking, V^O 



c -1 c - e -«S fej) 

Q 2 = SiS 2 e* 2 -=7 — c A2 =-sis 2 . 



Q 2 = , c A2 = -/(«, . 



Appendix B: Wave vectors of zigzag bilayer carbon tubes 

For AB-a stacking and V = the possible values of k^ v . are solutions of one of the equations 

3 cosQkMa) 3 sin(i K ( 2 )Af) 



or 



i + C ° 8 v 2;V sin(i K (DA) +S3 cos(i^)A) J " ° 
When V 7^ then the equation for k reads 



(1) sin(l K W(A + l)) (2) sin(^( 2 )(AT + l)) 
' * 3 sin(i K (DiV) * 3 sin(i«(2)A) 
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Here the function 



(E+vr-\0( K ,o\* 



describes the contribution of the second sheet of graphene to the eigenvector. 

For AB-/3 stacking and V — the possible values of K are solutions of the equation 

_ sin(.a) (j V + l)) \ / 2) /A sin(^) (j V + l)) \ 

sin( K (D7V) ) I + 3 \2j sin(«(2)7V) J 



' + 2S 3 ' COS (Jj .;„,.,,,, v.''' ! I 1 - 

(1) (2 

2 Sa S3 V^V 2 ^""1 2 y ' sin(K«iV)sin(K;<: 2 )iV) y 2 J V 2 77 2 



1 (i) (2) „.V« (2)> \ , 1 - cos(K«iV) cos^iV) ^ ^k (1)> \„._ f^ 2) \\ 1 



+ s r («» ^— J cos ( — ) - ; ;,, v ; v ..; ^ [ — ) si,. ( — ) j - - = o ibm 

When 7^0 then the equation for n is 



1 - cob(kWJV) cos(k( 2 )JV) . /k^ 1 ^ /«® . . . 

■sm — sin — =0 (B6) 



sin(/t( x )iV)sin(«;( 2 )Ar) V 2 / \ 2 
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